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We report electronic structure calculations for BaTi2Sb2 and discuss the results in relation to 
the observed superconductivity of this material when hole doped with Na. The Fermi surface shows 
several sheets. These include a nested nearly 2D cylinder. There are also two sheets, which are 
three dimensional, in spite of the layered crystal structure. A magnetic instability associated with 
Fermi surface nesting is found. A sign-changing s-wave state, different from the one in the Fe-based 
superconductors, is predicted within a scenario of spin-fluctuation mediated superconductivity. 



The discovery of high temperature superconduc- 
tivity in Fe-based compounds^ has lead to renewed 
interest in non-oxide superconductors. Recently, 
Doan and co-workers^ reported superconductivity in 
Bai_^Na^Ti2Sb20, 0.05<x<0.33 and up to 5.5 K. 
The structure of the parent compound, BaTi2Sb20, as 
determined by x-ray refinement,^ is shown in Fig. [I] 
This structure shares a number of common features with 
the Fe-based superconductors. Specifically, it is a layered 
structure containing square planes of transition element 
atoms, in this case Ti. Like the Fe-based superconduc- 
tors, there are two transition metal atoms per chemical 
unit cell. As may be seen, the two Ti in the BaTi2Sb20 
unit cell are related by a translation of [1/2,1/2] followed 
by a 90° rotation, and not by the glide of the Fe2As2 lay- 
ers in the Fe-based superconductors. Superconductivity 
was also very recently reported by Yajima and co-workers 
in undoped BaTi2Sb20 with a bulk of 0.8 K.^ 

The nearest neighbor Ti-Ti distance is 2.91 A, which is 
only somewhat larger than the nearest neighbor distance 
of 2.51 A in Ti metal, implying the possibility of impor- 
tant direct Ti-Ti bonding, as for the Fe lattice in Fe-based 
superconductors .1^ Also similar to the Fe-based supercon- 
ductors each Ti atom is coordinated by four pnictogens, 
although in the present case the geometry is that of ap- 
proximate square planes with Ti at the center, instead 
of approximate tetrahedra and furthermore, no Fe-based 
superconductor is known where the pnictogen is Sb. Fur- 
thermore, the Ti are also coordinated by O in the present 
compound. Each Ti has two O nearest neighbors. Thus 
as shown in the right panel of Fig. [T] if one considers only 
Ti-0 hopping, one obtains a bipartite lattice consisting of 
two independent sets of one dimensional ...Ti-O-Ti-0... 
chains, running respectively along x and y. Similarly if 
one considers only Ti-Sb hopping one also obtains one 
dimensional chains, ...Ti-Sb2-Ti-Sb2..., running along y 
and X where Sb2 denotes an Sb above the plane and an 
Sb below the plane. Such ID objects might therefore 
form the basis of a strongly nested electronic structure. 
However, as mentioned, the short Ti-Ti nearest neigh- 
bor distance suggests metal-metal bonding, which would 
couple the two sub-lattices. Furthermore the Sb atoms 
are arranged in linear chains along the c-axis with Sb-Sb 
distances of ~c/2=4.05 A, which is short enough that c- 




FIG. 1. Layered structure of BaTi2Sb2 (left) and as viewed 
along the c-axis with the Ba spheres suppressed (right). 

axis hopping could be important, especially assisted by 
the Ba. Returning to the metal-metal bonding, one sees 
that the bonds run along the diagonal directions in the 
two Ti atom unit cells, i.e. forming a square network ori- 
ented at 45° to the ...Ti-O-Ti-0... and ...Ti-Sb2-Ti-Sb2... 
chains. 

Experiments^ show that stoichiometric BaTi2Nb20 
has a transition at T5=54 K, similar to several other 
related compounds .I^M^ This transition has strong signa- 
tures in resistivity, susceptibility and specific heat. The 
microscopic mechanism for this transition has not been 
demonstrated, but it is described as either a charge den- 
sity wave (CDW) or a spin density wave (SDW). This 
is consistent with electronic structure calculations for 
the compound Na2Ti2Sb20, which show a nested Fermi 
surface. 

In any case, as Na is substituted for Ba in 
NaajBai_a^Ti2Sb2 0, Tq is systematically suppressed and 
superconductivity arises with Tc of up to 5.5 K.^ Impor- 
tantly, superconductivity co-exists with the phase corre- 
sponding to Ts including near the maximum Tc. Super- 
conductivity occurs for 0.05<x<0.33, which nominally 
amounts to average Ti valences between Ti^-^^^+ and 
rp^3.i65+ aggijijiing that Sb occurs as Sb^~. Titanates 
with valences between 4+ and 3+ are known to have 
electron-phonon superconductivity as in the case of elec- 
tron doped SrTi03,^but they can also be magnetic and 
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exhibit non-trivial correlation effects as in e.g. TiOCl 
and LaTiOs.^^ 

Here we report first principles studies of BaTi2Sb20 
in relation to magnetism and superconductivity. We find 
that like its sister compound Na2Ti2Sb20, the Fermi sur- 
face is nested. This nesting leads to peaks in the suscep- 
tibility that are incommensurate but near the X points. 
These peaks are strong enough to result in an actual mag- 
netic instability, which we confirm by direct calculations. 
This shows that is most likely an SDW instability. 
Therefore the superconductivity of doped BaTi2Sb20 is 
in proximity to an antiferromagnetic SDW phase. This 
leads to a specific prediction of the pairing state in a 
spin-fluctuation mediated scenario, in particular a sign 
changing 5-wave state, but not the same one as in the 
Fe-based superconductors.^^ 

The electronic structure was calculated using standard 
density functional theory based on the Perdew, Burke, 
Ernzerhof generalized gradient approximation.^^. For 
this purpose we used the general potential linearized aug- 
mented planewave method including local orbitals^^ as 
implemented in the WIEN2k code.^^ The LAPW sphere 
radii were 2.4 bohr for Ba and Sb, 2.15 bohr for Ti 
and 1.7 bohr for O. The calculations are based on the 
experimental crystal structure.^ This structure contains 
one internal coordinate, zsb=0.2514 corresponding to the 
Sb height. The value obtained from total energy mini- 
mization of zsb =0.2522 is very close to the experimen- 
tal value. This is in contrast to the large discrepancies 
found in the Fe-based superconductors in standard non- 
magnetic density functional calculations.'!^ We used well 
converged LAPW basis sets, including additional local 
orbitals, along with dense Brillouin zone samplings, as 
high as a 48x48x24 mesh for the primitive tetragonal cell. 

The calculated band structure is shown in Fig. |2]and 
the corresponding electronic density of states (DOS) is 
shown in Fig. |3] along with projections. These calcula- 
tions were done relativistically, including spin-orbit, but 
spin-orbit does not have a large effect. It does, however, 
affect the details at Ep- A comparison of the DOS with 
and without spin-orbit is given in Fig. [4j 

The value of the DOS at the Fermi energy (with spin- 
orbit) is N{Ef)=A.2^ eV~^ on a per formula unit ba- 
sis. This corresponds to a bare specific heat 7=10.0 
mJ mol~^K~^ and bare Pauli susceptibility Xo(q = 
0)=1. 37x10"^ emu/mol. The calculated Drude plasma 
energies are hVtp^xx='^-^^ eV and hVtp^zz='^-^^ eV for the 
in-plane and c-axis directions, respectively. For isotropic 
scattering this implies a resistivity anisotropy of ~1.6, 
meaning that from a transport point of view this is 
a very 3D material. This is reminiscent of what was 
found previously by Pickett for the related compound 
Na2Sb2Ti20.io 

The density of states shows that the O 2p derived 
bands lie between eV and -5 eV, while the Sb 5p 
bands are in the energy range -4 eV - -1 eV, all with 
respect to Bp- Therefore, these shells are nominally full. 
Also there are no Ba derived occupied valence bands. 



> 




FIG. 2. Calculated band structure of BaTi2Sb20 including 
spin orbit. The Fermi energy is denoted by the horizontal line 
at eV. 




FIG. 3. Electronic density of states and projections onto the 
LAPW spheres. 



Thus Sb is a trivalent anion in this compound and Ti is 
in its trivalent d} state. The density of states around Ep 
is derived primarily from Ti d states, hybridized with Sb 
p. 

The orbital character is illustrated in Fig. |5j which 
shows the band structure in a 1 eV range around Ep 
emphasizing the character of different (i-orbitals in a so- 
called "fat bands" scheme, where the bands are shown 
with symbols having a size that is a small value to make 
all bands visible plus an enhancement of the size propor- 
tional to the projection of given orbital character onto the 
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FIG. 4. Comparison of the density of states around Ep with 
and without spin orbit. 



LAPW sphere. As may be seen, there are three d orbitals 
- (^^2, dx2_y2 and dxy - that contribute significantly at 
Ef. There is strong dxy character in this mixture. 

The Fermi surface is shown in Fig. [6] and the depen- 
dence on band filhng is given in Fig. [7[ There are three 
sheets of Fermi surface. The first is a very two dimen- 
sional, square electron cylinder around the zone center. 
As seen in the fat bands plot, the band making up this 
cylinder contains a mixture of the three d orbitals that 
participate at Ep. There is a three dimensional complex 
shaped electron section around the T — Z line. This sec- 
tion has dz^ character. Finally, there is a large three di- 
mensional hole section around X. Like the square cylin- 
der electron section, this part of the Fermi surface has 
mixed character derived from the three active d-orbitals, 
though the dominant character is d^^ . Since the electron 
count is even, this section compensates the two electron 
sections. This basic structure is preserved as the Fermi 
level is lowered within the range corresponding to su- 
perconductivity in Naa^Bai_ajTi2Sb20. The hole section 
grows, while the two electron sections shrink. Impor- 
tantly, the square cylindrical shape of the 2D section is 
preserved. 

An important aspect of the Fermi surface is that the 
square section, as well as the edges of the X centered elec- 
tron section are nested along the [1,1] direction. The 2D 
nesting vector is approximately (0.24,0.24) (27r/a). [1,1] 
is the orientation of the Ti square lattice. In this lat- 
tice, the dxy orbitals are along the Ti-Ti bond directions, 
and may be expected to give nested bands. Actually, 
the nested electron section has strong dxy^ but also as 
mentioned involves the other two active d orbitals. In 
any case, qualitatively the nesting found should lead to 
two ridges in the bare susceptibility xo running along 
the diagonals in the zone. The intersection near the X 
points would give a prominent peaks in Xo(q) near the 
X points. This in fact is the case. Fig. [S] shows the 
bare Lindhard xo calculated with the neglect of the ma- 
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FIG. 5. Fat bands plot for different d orbital characters near 
Ef- The coordinate system is such that d^^ is directed along 
the c-axis, and x and y are along the a- axis and h axis, such 
that dxz is towards the nearest Sb and y is in the direction of 
the O neighbors. 



trix element approximation, similar to calculations pre- 
sented for the Fe-based superconductors.'^ As seen there 
are diagonal ridges in the real part of xo^ with a promi- 
nent peak at X. This peak will be further increased by 
the RPA enhancement, x(q)=Xo(q)/(l-^(q)Xo(q)) (here 
/(q) is written as a Stoner term, but it in correlated ma- 
terials it takes the form of the Hubbard parameter U). 

The real part of x governs magnetic ordering as well 
as providing the pairing interaction in spin-fluctuation 
induced superconductivity.fi^^ We discuss it further be- 
low. The imaginary part of xo^ like the real part, shows 
ridges with intersections near X, but the structure is 
more complex. This structure in both the real and imag- 
inary parts of X, ^-e. intersecting ridges with peaks at the 
intersections, is a consequence of the Fermi surface struc- 
ture, particularly the square cylinder shape. This Fermi 
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FIG. 6. Calculated Fermi surfaces of BaTi2Sb2 spin orbit. 
The shading is by velocity, with blue representing low ve- 
locity. The approximately square cylindrical sections around 
the zone corner (M,A) are electron sections as is the complex 
shaped section around F and Z. The other section around R 
is a hole section. 
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FIG. 7. c-axis view of the Fermi surface as in Fig. [6] but 
with shifted Fermi energies. The corresponding band fillings 
are -0.38 e (-0.10 eV), -0.23 e (-0.05 eV), e (0 eV), and 0.20 
e (0.05 eV). 



surface amounts to the intersection of two ID sections. 
The result is similar to the antiferromagnetic peak in the 
triplet superconductor Sr2Ru04p^^ but different from 
the Fe-based superconductors where the nesting arises 
from matching of electron and hole Fermi surface sec- 
tions. 

We did direct self-consistent calculations to look for a 
magnetic ground state. We find a magnetic instability 
at the X point. These calculations were done includ- 
ing spin-orbit with the moment directed along the c-axis. 
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FIG. 8. Calculated c-axis projections of the real and imagi- 
nary parts of xo (see text) with arbitrary units. Light shading 
denotes higher values. 
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FIG. 9. Magnetic structure of the X-point ordering showing 
the two degenerate cases in the left and right panels. The 
doubling of the unit cell is shown. Red and blue denote up 
and down spins, respectively. Note the alternation of ferro- 
magnetic and antiferromagnetic bonds in this double stripe 
structure. 



We find an instability, with small Ti moments inside the 
LAPW spheres of 0.2 jub- This is much smaller than 
the ^ 2 jiB characteristic of the Fe-based superconduc- 
tors in similar calculations using the experimental crystal 
structure. We also checked for a nearest neighbor an- 
tiferromagnetic state, which would correspond to a zone 
center instability. However, we did not find such a state. 

The magnetic structure is shown in Fig. [9j This a 
so-called double stripe structure. It lowers the symmetry 
from tetragonal to orthorhombic and contains both ferro- 
magnetic and antiferromagnetic Ti-Ti bonds, as shown. 
This may be expected to result in a spin-dimerization, 
with coupling to the lattice, perhaps consistent with neu- 
tron scattering results showing lattice anomalies at Tg 
in Na2Ti2Sb20.^ The actual magnetic structure may be 
incommensurate, as the moments are small and as men- 
tioned, the intersection of the ridges in x(q) is off the X 
point. Another indication that the actual ground state 
may be a more complex incommensurate magnetic struc- 
ture is that in our calculations without imposed symme- 
try for the X point instability we find slightly different 
moments on the two Ti sites of each spin (the magnetic 
unit cell contains four Ti atoms, two spin-up and two 
spin-down) . 
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In any case, we find a significant coupling of the mag- 
netic order to the electrons at the Fermi energy. In partic- 
ular, even though the moments are small, the DOS near 
Ef is substantially reduced in the ordered state. We find 
N{Ef)=3A2 eV~^ per formula unit for the X point or- 
dered state as compared to N{Ef)=4:.26 eV~^ without 
magnetic order. As in the Fe-based superconductors,'^ 
and in contrast to cuprates, the antiferromagnetic state 
near superconductivity is metallic, both in experimenlP 
and in the present density functional calculations. 

One of the reasons why unconventional superconduc- 
tivity attracts so much attention is that it is a rare phe- 
nomenon. In particular, the vast majority of known 
superconductors are conventional electron-phonon su- 
perconductors, and in particular conventional electron- 
phonon superconductivity is known in titanates. Unlike 
the Fe-based superconductors, the maximum Tc=5.5 K 
observed in BaTi2Sb2 is not so high as to preclude con- 
ventional superconductivity. However, the close prox- 
imity to an antiferromagnetic state provides an indica- 
tion that this may not be the case. In particular, spin- 
fiuctuations provide a repulsive interaction in a singlet 
channel, opposite to the attractive electron-phonon in- 
teraction. Therefore the net superconducting interac- 
tion will be A = Xep — A^/, where ep and sf denote 
electron-phonon and spin-fiuctuation contributions, re- 
spectively. As such the wide coexistence, of magnetic 
and superconducting order in the composition depen- 
dent phase diagram argues against (but does not fully 
exclude) conventional superconductivity. An argument 
against an unconventional state would be that scattering 
due to disorder in the alloy should suppress an unconven- 
tional state. However, the disorder in the present case is 
in the Ba layer, and as mentioned Ba does not contribute 
significantly to the electronic structure near Ef^ and fur- 
thermore, as mentioned the bands are rather fiat, which 
means that the superconducting coherence length will be 
relatively short, again reducing the effect of disorder. 

As such, we discuss the possibility of a spin-fiuctuation 
mediated unconventional superconducting state. In such 
a scenario, the pairing interaction is related to the real 
part of the actual RPA enhanced x(q) is re pulsive 
for a singlet state and attractive for a triplet state.^i^*^^^ 
As mentioned, in BaTi2Sb20 the susceptibility is peaked 
near X. Identification of the superconducting state 
amounts to matching the pairing interaction with the 
Fermi surface structure, so that opposite sign order pa- 
rameters occur on parts connected by peaks in the repul- 
sive (singlet case) interaction or parts connected by peaks 
in the attractive (triplet case) interaction have the same 
sign order parameter, while maintaining consistency with 
the overall triplet or singlet parity. We find that the best 
match is obtained for a singlet state as depicted in Fig. 
[Tol In this state the two sections of Fermi surface with 
similar orbital character, namely the sqaure section and 
the large section around X, have opposite sign as ex- 
pected in such a scenario. 

As seen, this is a sign-changing 5-wave state, although 




FIG. 10. Potential sign changing s-wave superconducting 
state, shown on top of the Fermi surface. Note that in this 
figure, r is at the center, as opposed to the corners of the 
plot. 



not the same one as in the Fe-based superconductors. It 
consists of a state where the two electron sections of the 
Fermi surface have opposite sign order parameter to the 
large hole section. The motivation for this state is that it 
has sign changes between all pairs of Fermi surface sec- 
tions that are separated by (0,7r/a) and (7r/a,0), where 
the susceptibility is peaked. Similar to the Fe-based su- 
perconductors, the SDW and this superconducting state 
are competing instabilities of the Fermi surface. How- 
ever, the SDW antiferromagnetism found for BaTi2Sb20 
is much weaker than that in the Fe-based superconduc- 
tors. It is unclear what role correlation effects beyond 
standard density functional calculations play in these ti- 
tanates. The fact that the transition at Tg is a metal- 
metal transition in spite of the integral band filling sug- 
gests that they are not so strong as in e.g. the cuprates. 
This of course is an expected result considering that the 
d orbitals in Ti^+ ions are more extended than those in 
Cu^+. One possible effect of Coulomb correlations will 
be to enhance the small moments found here. It will 
be of interest to perform neutron diffraction experiments 
both to confirm whether the transition is in fact an SDW 
and to quantify the magnitude of the moments. In this 
regard, Subedi has recently presented electron-phonon 
calculations^^ that are consistent with electron-phonon 
superconductivity and so either scenario is possible here. 

Returning to Na2Ti2Sb20, that material, which has 
been studied more extensively than BaTi2Sb2 0, has not 
been reported to show superconductivity, but does have 
a transition, likely of spin density wave or charge den- 
sity wave character and a similar phenomenology to 
BaTi2Sb20. The transition in Na2Ti2Sb20 is at higher 
temperature, Ts~114 K - 120 K.^ ^ From a structural 
point of view, Na2Ti2Sb20 has similar Ti2Sb2 layers 
to BaTi2Sb20, but these are stacked differently along 
the c-axis so that the Sb-Sb linear chains in the Ba com- 
pound are broken up yielding a body centered tetragonal 
structure. The band structure, which as mentioned has 
similarities to the present compound, was investigated 
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by Pickett using density functional theory, and in a tight 
binding framework by de Biani and co-workers. The 
arsenide, Na2Ti2As2 has a transition at Ts~320 K.^ 
Other related compounds may be (SrF)2Ti2Sb20, which 
has an apparently similar phase transition at T^^IOS K 
and (SmO)2Ti2Sb20 with T,-230 KPBaTi2As20 also 
has been reported and has similar transition at T5~200 
K. Considering the superconductivity of BaTi2Sb20, it 
will be of great interest to study these compounds in 
more detail and especially their properties if the SDW 
can be suppressed, e.g. by pressure, doping or chemical 
substitutions. It will be particularly interesting to study 
the pressure dependence in BaTi2Sb20 and related ma- 
terials. 

In conclusion, we find that BaTi2Sb2 has an SDW in- 
stability associated with Fermi surface nesting. The close 
proximity of the superconducting phase that emerges 



with Na doping to this antiferromagnetic phase suggests 
the possibility of unconventional spin-fiuctuation medi- 
ated superconductivity. We find that the susceptibility 
is peaked near the X point. Matching this structure of 
the susceptibility with the Fermi surface suggests the pos- 
sibility of a sign-changing 5-wave superconducting state 
of a different nature than the Fe-based superconductors. 
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